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BOSE-EINSTEIN CONDENSATION OF FERMION PAIRS 
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Abstract. We consider the low density limit of a Fermi gas in the BCS approxi- 
mation. We show that if the interaction potential allows for a two-particle bound 
state, the system at zero temperature is well approximated by the Gross-Pitaevskii 
functional, describing a Bose-Einstein condensate of fermion pairs. 



1. Introduction and Main Results 

1.1. Introduction. The bosonic behavior of pairs of fermions is a topic that has 
been investigated in condensed matter physics for more than half a century. It plays 
a crucial role in the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity [2], 
and is used to explain the superfluid behavior of He^ and cold gases of fermionic atoms, 
for instance. As long as the pair of fermions is tightly bound, it is not surprising that 
it effectively behaves like a boson, and hence can form a Bose-Einstein condensate 
(BEC). In BCS theory the pairing mechanism is also important under very weak 
attraction where the separation of the paired particles can be much larger than the 
average particle spacing, however. 

It was realized in the eighties [141 [15] that BCS theory actually applies both to the 
case of BECs of tightly bound fermions and to cases where the pairing mechanism is 
very weak. The regime in-between is called the BEC-BCS crossover regime |18[ H]. This 
crossover is currently a hot topic in condensed matter physics, and recent experiments 
on cold atomic gases have been able to probe large parts of this regime. We refer to 
[3] for a recent review. 

From the mathematical physics point of view, the pairing mechanism in fermionic 
systems is poorly understood, and there are no rigorous results starting from first 
principle, i.e., with an appropriate many-body Hamiltonian. In this work, we shall 
assume the BCS approximation to be correct, and investigate some of its consequences. 
This paper can be viewed as a sequel to the recent work |6l where the emergence of 
Cinzburg-Landau (CL) theory [8] from BCS theory was studied. Close to the critical 
temperature, GL arises as an effective theory on the macroscopic scale, describing the 
variations in the density of fermion pairs. For this it is not necessary to form actual 
bound states between the fermions, a very weak attraction is sufficient for pairing. 

In this paper we are interested in the low density limit in the case where the inter- 
particle interaction does allow for two-particle bound states. We consider the system 
at zero temperature and show that the macroscopic variations in the pair density are, 
to leading order, correctly described by the Gross-Pitaevskii (CP) functional [9l [T7]. 
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describing a BEC of fermion pairs with an effective repulsive interaction. We allow 
for a large class of possible interactions among the particles. Essentially our sole as- 
sumption will be the existence of two-particle bound states. The existence of such 
bound states is crucial here. In their absence one obtains an ideal Fermi gas in the 
low density limit, up to exponentially small corrections [5l [121 113) . 

The proof of our results uses the same tools as the proof of the main theorem in 
[B]. Parts of it are simpler, in fact, since we work at zero temperature here. We shall 
demonstrate that the semiclassical estimates in [6] extend to the zero temperature 
case. 

1.2. The BCS Functional. We consider a macroscopic sample of a system of spin 1/2 
fermions at zero temperature. For simplicity, we restrict our attention to three spatial 
dimensions, but our analysis applies to any dimension 1 < d < 3. The interaction 
among the fermions is described by a local two-body potential V. In addition, the 
particles are subject to external electric and/or magnetic fields. Neutral atoms would 
not couple to these fields, of course, but there can be other forces, e.g., arising from 
rotation, with a similar mathematical description. In BCS theory the state of the 
system is described in terms of a 2 x 2 operator valued matrix 

- /-) M 
a 1 — 7 / 

satisfying < F < 1 as an operator on © L^(M^). The bar denotes complex 



conjugation, i.e., a has the integral kernel a{x, y). The fact that F is hermitian implies 
that 7 is hermitian and a is symmetric, i.e., 7(3;, y) = ^{y,x) and a{x,y) = a{y,x). 
Moreover, since F^ < F, we have < 7 < 1 and < aa < 7(1 — 7). 

We are interested in the effect of weak and slowly varying external fields. Hence 
we introduce a small parameter h > and write the external magnetic and electric 
potentials as hA{hx) and h'^W{hx), respectively. In order to avoid having to introduce 
boundary conditions, we assume that the system is infinite and periodic with period 
h^^, in all three directions. In particular, A and W are periodic, and we assume that 
the state F is periodic. Our goal is to calculate the ground state energy per unit 
volume, and the corresponding BCS minimizer. 

We find it convenient to do a rescaling and use macroscopic variables instead of the 
microscopic ones. The rescaled BCS functional has the form 



£^^^{T) := Tr ( {-ihV + hA{x)y - + h^W{x)j 7 

V{h'\x-y))\aix,y)\Uxdy (1.2) 



Cxi 

where C denotes the unit cube [0, 1]^, and Tr stands for the trace per unit volume. 
Explicitly, if x denotes the characteristic function of C, and B is a periodic operator 
with x^X trace class, Tr B equals the usual trace of x^X- The location of the cube is 
obviously of no importance. Using the Floquet decomposition [191 Sect. XIII. 16], it is 
not difficult to see that the trace per unit volume has the usual properties of a trace 
like cyclicity, for instance, and standard inequalities like Holder's inequality hold. This 
is discussed in detail in [6l Sect. 3]. 
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In ()1.2p we choose units such that the particle mass equals 1. The particles have 
spin 1/2, which adds an extra factor 2 to the energy. The chemical potential is denoted 
by /i/2, for convenience, and the external electric potential is really W 12. 

For heuristic arguments explaining the derivation of the BCS functional ()1.2p . we 
refer to [101 Appendix A]. The BCS state of the system is a minimizer of this functional 
over all admissible states F, i.e., periodic F of the form (11. ip satisfying < F < 1. 

We make the following assumptions on the potentials A and W in ()1.2p . Our results 
presumably hold under slightly weaker regularity assumptions on W and A, but to 
keep things simple we shall not aim for the weakest possible conditions. 

Assumption 1. We assume both W and A to be periodic with period 1. We further 
assume that W{p) and \A{p)\{l + \p\) are summable, with W{p) and A{p) denoting 
the Fourier coefficients of W and A, respectively. In particular, W G C''(]R^) and 
AeCU 



The interaction potential will be assumed to satisfy the following properties. 

Assumption 2. The interaction potential V is assumed to be real-valued and reflection- 
symmetric, i.e., V{x) = V{—x), with V £ L^/^(R^). Moreover, the Schrodinger oper- 
ator — + V{x) has a negative energy bound state. 

The L3/2 assumption on V guarantees relative form-boundedness with respect to 
the Laplacian. The ground state energy of — + V{x) will be denoted by — .Bb < 0, 
and its ground state wave function by ao- It is unique up to a phase factor. We find 
it convenient to normalize oq such that 

dq 
(2^3 

with ao(g) = (27r)~'^/-^ J^^ a{x)e~^'^'^dx denoting the Fourier transform. 

In the following, we are interested in the case ^ = —E\y + h'^Sfx, which corresponds 
to the low density limit. We find it convenient to absorb the constant h^S/i into the 
potential W, i.e., we set /i = — -Eb and write W{x) instead of W{x) — 6fj,. 

1.3. The GP Functional. Let -0 e H^^^{R^), the periodic functions in Hl^^{R^). For 
g > 0, the GP functional is defined as 

= / \l\{-iV + 2A{x))ij{x)\'^ + W{x)\^p{x)\'^ +g\i;{x)\^]dx. (1.4) 



\Mq)\'7^ = i, (1-3) 



c 

The coefficient 2 in front of the vector potential A is due to the fact that ^ describes 
pairs of particles, and the charge of a pair is twice the particle charge. The factor 4 in 
front of the kinetic energy is twice the mass of a fermion pair. The coefficient g will 
be calculated below from BCS theory. 

We denote the ground state energy of the GP functional by 

E'^^'ig) = inf {£{^P) : V G H^.^iR')} . (1.5) 

It is not difficult to show that under our assumptions on A and W, there exists a corre- 
sponding minimizer, which satisfies a second order differential equation known as the 
GP equation. Note that there is no normalization constraint on ip, the normalization 
is determined by the chemical potential which is contained in W{x). 
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aoiq)\Hq' + E^)j^. (1.7) 



1.4. Main Results. We define the energy E^''-^^{fi) as the infimum of S^*^^ over all 
admissible F, i.e., 

EBCS(^)^i^f ^BCS(r). (1.6) 

Recall that a state F is admissible if < F < 1 and F is periodic, i.e., it commutes 
with translations by 1 in the three coordinate directions. 

Recall also that oq denote the ground state of — + V{x), normalized as in (II. 3p . 

THEOREM 1. Let 

dq 
(27) 

Under AssumptionsUl and\E above, we have, for small h, 

E''''^{-E^) = h{E''^{g)+e) , (1.8) 

with e satisfying the bounds const /i > e > — const /i"*^/^. Moreover, ifV is an approxi- 
mate minimizer of £^^^ at fj, = —E]^, in the sense that £^^'^^(F) < h{E^^{g) + e) for 
some small e > 0, then the corresponding a can be decomposed as 

a = ^{7p{x)ao{-ih\/) + aQ{-ihV)Tp{x)) + a (1.9) 

with £'^^{^) < E'-^^{g) + e + const /i^^ and 

\a{x,y)\^ dxdy < const h^/^ . (1.10) 



To appreciate the bound (jl.lOp . note that the square of the L^(C x R^) norm of the 
first term on the right side of (II. 9p is of the order h~^, and hence is much larger than 
the one of a. To leading order in h, the pair wave function a{x, y) is thus given by 

ip{x)+'ip{y) 1 
with ip a minimizer of the GP functional (jl.4p . This agrees with 

to leading order in h, hence tp describes the center of mass motion of pairs of fermion 
with are bound in the ground state of — /i^V^ + V{x/h). 

Note that g in (jl.7p is strictly positive, even for purely attractive interaction poten- 
tials V . In the limit of a point interaction [U Sect. LI] with scattering length > 
we have 

^El" , ^ 1 



So(g) = o , ^ and ^Jfe = — , (1.13) 

and hence g = liras. Since the mass of the fermion pairs is 2, this corresponds to a 
scattering length of 2 Ug for the pair scattering [20] . The factor 2 is an artifact of the 
BCS approximation; an investigation of the actual four-body problem with pseudo- 
potentials predicts a scattering length 0.6 |16j . 

By varying the external potential W, our bounds (II. Sp on the ground state energy 
can be used to obtain bounds on the particle density as well. In particular, the number 
of particles per unit volume, A^, can be calculated by replacing W{x) by W(x) + 6fi 
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and taking the derivative of the energy with respect to ^h?5^ at dfi = 0. To leading 
order in h, the result is that 

N = l [ \^p^^{x)fdx (1.14) 
" Jc 

where ip^^ is a minimizer of the GP functional (jl.4p . The average particle density, in 
microscopic variables, is p = h^N = 2h? f and is thus of order h"^. Hence our 

scaling limit corresponds indeed to low density. 

In the translation invariant case, where A(x) = and W{x) = —Sfi < is constant, 
the GP minimizer is given by |'0*^^(x)p = 6p/{2g), and E^^{g) = —6p?/{4,g). In 
particular, the ground state energy per particle, which is equal to E^'^^{fi)/N + ^fi, 
is given by 

— + ^gp + higher order in p (1-15) 

for small density p. 

1.5. Outline of the paper. In the following Section [2] we shall state our main semi- 
classical estimates. These are a crucial input to obtain the bounds in Theorem [H 
They are an extension to zero temperature of the analogous expressions at positive 
temperature obtained in pi Sect. 2]. An upper bound on E^^^ will be derived in 
Section [3l using the variational principle. Finally, Section U] contains the lower bound. 
In this final section also the structure of approximate minimizers will be investigated. 
This leads to a definition of the order parameter ip. Our proof follows closely the proof 
of the main theorem in [6], but is partly simpler due to the fact that we work at zero 
temperature. 

Throughout the proofs, C will denote various different constants. We will some- 
times be sloppy and use C also for expressions that depend only on some fixed, h- 
independent, quantities like E^ or ||VF||oo; for instance. 



2. Semiclassical Estimates 

This section contains the semiclassical estimates needed in the proof of Theorem [TJ 
Let tp he a periodic function in Hf^^{M.^). Pick a reflection-symmetric and real- valued 
function t, with the property that 

d'^t e l6(r3) (2.1) 

and 

(2.2) 

1 + 

for all multi-indices 7 G {0, 1 . . . ,4}^. We shall later choose t{q) = 2{q^ + £'b)So(Q'), 
but the results of this section are valid for general functions t satisfying (j2.ip and (j2.2p . 
Let A denote the periodic operator 

A = (V'(x)t(-iW) + t(-iW)V'(x)) , (2.3) 
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and let 

_ r {-ihV + hA{x) f - H + h^W{x) A 

^ -{ihV + hA{x)f + n-h'^W{x) 

(2.4) 

on L^(M^) C^, with A and W satisfying Assumption [H We shall also assume that 
fj, < 0. In the following, we will investigate the trace per unit volume of the negative 
part of H/^. Specifically, we are interested in the effect of the off-diagonal term A in 
H/\, in the semiclassical regime of small h. 

THEOREM 2. Let [s]_ = ^ {\s\ — s) denote the negative part. For /i < 0, the 
diagonal entries of the 2x2 matrix-valued operator [H/\]^ — [Hq]^ are locally trace 
class, and the sum of their traces per unit volume (which will be denoted by Ttq ) equals 



TtoHHa]- - [Ho]-) = -h-^E^ - hE2 + 0{h') (llV'll^i(c) + MmiO^ 

+ 0{h')[\m%r^c) + MW)) ' (2.5) 



where 



2 7r3 r - ^ (27r)- 



and 



3 

^2 = V (a,V|9fcV> / t{q) [djdkt]{q) 



dq 



2^"' '"VV(g2-/,)2(2vr)3 

The error terms in i2. 5]) of order h'^ and h^ depend on t only via bounds on the 
expressions 112. 1\) and 112. 2\) . and are uniform in /j, for /j, < bounded away from zero. 

Here, we use the short-hand notation HV'llp for the norm on LP(C). Likewise, ( • | • ) 
denotes the inner product on L'^{C). 

In general, the operator [H/s\- — [Hq\- is not trace class under our assumptions on 
t and ■0- Hence the trace in (|2.5p has to be suitably understood as the sum of the 
traces of the diagonal entries. This issue is further discussed in the next section. 

The proof of Theorem [2] is very similar to the proof of Theorem 2 in [6j . In the 
following, we shall limit ourselves to explaining the main differences. 

Sketch of proof. Since /i < and W{x) and '4){x)t{—ihV) are bounded, both //a and 
Hq have, for small enough /i, a gap around in the spectrum. Hence the projector 
onto the negative spectral subspace can be written via a contour integral as 

e{-HA) = ^. [ —^dz, (2.8) 
Z-Kl J ^ Z — ilA 

where 9{t) = 1 for t > 1 and otherwise, and i is the contour {r — i, r G (— oo,0]} U 
{ir, r G [—1, 1]}U{— r-|-i, r G [0,oo)}. The integral has to be understood as a suitable 
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weak limit over finite contours. Similarly, one obtains that 



1 

27ri 



1 



1 



z — z — 



dz . 



(2.9) 



The remaining analysis proceeds as in ^ Sect. 8] (compare with Eq. (8.11) there), 
and we shall not repeat it here. In [6], the factor z on the right side of (|2.9p is replaced 
by —f3~^ ln(l + e~^^) and the contour is around the whole real axis. For /3 — )• oo, this 
reduces to ()2.9p . given the gap in the spectrum around 0. □ 

Our second semiclassical estimate concerns the upper off-diagonal term of the pro- 
jection onto the negative spectral subspace of H/\, 9{—H/\), which we denote by a a.. 
We are interested in its norm. In general, we define the norm of a periodic 
operator O by 



\0\ 



Tr 



Ot (1 - /i2v2) O 



(2.10) 



In other words, ||0||^i = ||0||2 + /i^||VO||2. Note that this definition is not symmetric, 
i.e., ||0||hi / llO^II/fi in general. 

THEOREM 3. Let ip{q) = \t{q)/{q^ - fi). Under the same assumptions as in 
Theorem\^ we have 



{^{xM-ihV) + ^{-ihV)^{x))\\^, < Ch^l^ (UWhhc) + UfnHc) 



(2.11) 



The proof follows again along the same lines as the proof of the corresponding 
Theorem 3 in [6], and we shall only sketch the differences. 



Sketch of proof. With the aid of ()2.8p we can write 



OA 



1 

27^ 



z-H^ 



dz . 



(2.12) 



12 



where stands for the ij element of an operator- valued matrix, and where the 



integral has to be suitably understood as a weak limit, similarly to (j2.9p . Alternatively, 
one could integrate over —I, since the identity operator has vanishing off-diagonal 
terms. 

Using the resolvent identity and the definitions of A and ip we find that 

h ^ 
OA = - {ip{x)ip{-ihV) + ip{-ihV)ip{x)) + ^ r/j , 

i=i 

where 



(2.13) 



■ni 



- A 

Ani 
1 

2TTi 



1 



Z — kr 



z — k, 



and 



'0 

1 

2^i 



{k - ko) 



t 



kt 



+ 



t 



1 



k^ z + ko 



dz 



' A + A ' 



-A 



z — k 
1 



z + h 



■(ko-k) 



1 



■At- 



z — k z + k z — k 



-A 



z + k 
dz . 



dz 



(2.14) 
(2.15) 

(2.16) 



22 
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Here, t is short for the operator t{—ihV), k = {—ihV + hA{x))'^ — A* + h?W{x) and 

Proceeding as in ^ Section 9] one sees that 

\\vi\\m<Ch^^''\mHHc), (2.17) 
\\m\\m<Ch'^^\mm(c) (2.18) 

and 

\\V3\\m<Ch'/^^pfH^^e)■ (2.19) 
In the terms investigated in [6j, there is an additional factor (l+e^^)~^ in the integrand, 
and the contour contains the whole real axis. Similarly as in the proof of Theorem [21 
this reduces to our case as /3 — >• oo. □ 

3. Proof of Theorem [1} Upper Bound 

Recall that oq denotes the unique ground state of — + V{x), normalized as in 
()1.3p . It satisfies ao{x) = ao{—x), and we can take it to be real. In the following, we 
let t denote the Fourier transform of 2(— + E]^)ao = —2^00, i.e., 

t{q) = -2(2^)-3/2 / y(x)ao(x)e-^«-"dx = 2{q^ + ^b)3o(g) . (3.1) 



It satisfies all the assumptions in the previous section. In particular, (j2.ip and (|2.2p 
hold for all 7 G Nq. This can be shown, for instance, in the same way as in [6j Sect. 4]. 
The method there also implies that y^|F(x)|Qo(x)e'^l^l G L^(R'^) for k < E^"^ , and 
that J^-i{\x'^'S/ ao{x)\^ + \x'^ ao{x)\'^)dx < 00 for all 7 G Ng. Some of these properties 
will be used later on. 
As a trial state, we use 

Fa = f !^ )= ei-HA) (3.2) 

V OA 1 - 7A / 



where H/^ is given in (j2.4p with A as in (j2.3p and fj. = —E^,. For t, we choose ()3.ip . 
which is refiection symmetric and can be taken to be real. 
We have 

[i^o]- - [Ha]- = HaTa - HoTo = ( '.^^ ) (3.3) 

V ^7A + kaA K7A + AoA / 

where k denotes the upper left entry of Ha (and -f^o). From ()1.2p and p.3p we conclude 
that 

fBCS(r^)^_ l^ro [[Ha]- - [Ho]-) 

-h-' i _^n^)||(V'(x)+V(2/))ao(^)|' ^ 



+ 



(2vr) 

f ^r/x-ys ^{x)+^{y) 

/cxMa"^^^^ 2/.^(2vr)3/^ «o(^) " aA(x, ,) 



2 

dx dy , (3.4) 



where Trg stands for the sum of the traces per unit volume of the diagonal entries of 
the 2x2 matrix- valued operator. In general, the operator [i?o]- — [-^^a]- is not trace 
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class if A is not, as can be seen from ()3.3p . In the evaluation of iS^'~'^(rA) only the 
diagonal terms of (j3.3p enter, however. 

The first term on the right side of (13. 4p was calculated in Theorem [2] above. Note 
that, for our choice of t, the integral in the second term in ()2.7p is equal to 4, and the 
integral in the third term is 16 g. 

As in [6', Sect. 5] we can rewrite the second term on the right side of ()3.4I) as 

[ V{h-\x-y))\^{i;{x)+i;{y))ao(,h-Hx-y))\^ dxdy 

Using the Taylor expansion 

2t{q) + t{q - hp) + t{q + hp) 



At{q) + h^ [{p ■ V)h] (q) + ^ J'^ [ip . V)H] (q + shp)il - \s\)^ ds 



(3.6) 



we see that (|3.5p equals 



4^dq+^ Y.{dMd,^) f t{q)[d,d,t]{q)^^dq + 0{h^)A^-^) 
J^z q^ + Ey, 16 J^3 q^+Ey, 



Ah 



where the error term is bounded by 



Ch^i;\\l, / \Vix)\\aoix)\^\x\Ux. (3.8) 



1^2 

Note that the first term in p.7p cancels the contribution of h^^Ei in (12. 6p to the trace 

iT:([i/A]--[i^o]-). 

It remains to investigate the last term in (j3.4p . Since V is relatively bounded with 
respect to the Laplacian, we can bound the term by an appropriate norm. Recall 
the definition of the norm of a periodic operator in (|2.10p . For general periodic 
operators O, we have the bound 

/ V{^)\0{x,y)\^ dxdy < {I -V^)-y^V{-)(l -V^)'^'^ \\0\\]j, . (3.9) 

The operator of relevance here is given by 

O = OA - I (V'(a^)ao(-i/iV) + ao(-i/iV)V'(x)) . (3.10) 

Note that, for our choice of t, we have ip{q) = \t{q)/{q^ + -Eb) = ao{q)- Hence 
Theorem [3] implies that the norm of ([330]) is bounded by Ch^/'^{\\il;\\H^ + \\i^\?Hi)- 
For ^jj, we shall take a minimizer of the GP functional (II. 4p . Under Assumption [T] 
on W and A, it is easily seen to be in H^. Collecting all the terms, we see that for 
this choice of "0 we have 

for small h. This completes the proof of the upper bound. 
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4. Proof of Theorem [TJ Lower Bound 

Our proof of the lower bound on -Eb) iii Theorem [1] consists of two main 

parts. The goal of this first part is to show the following. 

Proposition 1. Let T be a state satisfying E^^^iT) < 0, and let a denote its off- 
diagonal element. Then there exists a periodic function ip, with H^{C) norm bounded 
independently of h, such that 

a = |(V'(x)ao(-iW) + aQ{-ihV)^{x)) + i (4.1) 

with W^Whi < 0{h^/^) for small h. 

Recall the definition (|2.10p for the norm of a periodic operator. The bound 
IICIIhi < 0{h^/'^) has to be compared with the norm of the first part of (14. ip . 
which is 0(/i-^/2) (fo^. fl^ed / 0.) 

Proof. Let K"^'^ denote the operator 

kA,w ^ ^_-f^^ ^ hA{x)f + + h^W{x) , (4.2) 
and A"^'^ = inf spec X"^'^ > -E'b ~ 0{h?). For all admissible states L, we have 

< aa < 7(1 -7) < 7 (4.3) 

and hence 

Tr K^''^7 > TV K^'^ [aa + 7^) > TV K^^^ aa + A^'^TV 7^ 

> TV K'^'^aa + A^'^Tr {aaf . (4.4) 

In particular, 

^BCS(r) > ^A,W^, („^)2 + ^^^( . y-^\KA,W + _ . (4 5) 

Here K^'^ acts on the x variable of a{x,y), and (• | •) denotes the standard inner 
product on L^(M'^). 

By definition, the operator K^'^ + V{h~^{ ■ — y)) on L^(M'^) has a unique ground 
state h~^^'^ao{h~^{x — y)), with ground state energy zero, and a gap above. To utilize 
this fact, it will be convenient to replace K"^'^ by K'^'^ in (j4.5p . We claim that 

K^^w + y^f^-i^, -y))>\ (^°'° + Vih-\. - y))) - {\\W\U + PIlL) • (4.6) 
This follows immediately from the fact that 
K^,w ^ y^^-i^, -y)) = \ {K''' + V{h-\- - y))) 

+ ^ (/C2A,0 ^ . _ y))) _ h^Aix)^ + h^Wix) (4.7) 

and that K"^^'^ + V{h~^{ ■ — y)) > 0, by the diamagnetic inequality. 

For any state F with iS^*-'^(F) < we conclude from (j4.5p and (j4.6p that 

A^'^TV {aaf + ^ ^(a( • , y)\K'^^ + V{h-\ - - y))\a{ ■ , y)) dy < Ch^\\a\\l . (4.8) 

To show that this inequality implies (|4.ip . we shall proceed as in [6l Sect. 6]. 
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Define ip to be the periodic function 

ipiy) = i/n, / ao{h'\x - y))a{x,y)dx . (4.9) 

If we write 

<^{x,y) = j^;^^^^j^^ao{h~\x - y))ij{y) + ^o{x,y) (4.10) 

the gap in the spectrum of K^'^ + V{h~^{ ■ — y)) above zero, together with (j4.8p and 
the normahzation ()1.3p . yields the bound H^olb < O(^)||o||2- We can also symmetrize 
and write 

a{x,y) = '^2(2Tr)^th^ ao{h-\x - y)) + ^(x, y) , (4.11) 

again with ||^||2 < 0(/i)||a||2. In order to complete the proof of (j4.ip . we need to show 
that 1 1 "01 1 HI is bounded independently of h, and that the norm of ^ is bounded by 
0(/iV2). 

An application of Schwarz's inequality yields 



Pix)\^ dx < h\\ag < \i;ix)\^dx + h\Ml . (4.12) 
Since H^olb < 0(/i)||a||2, this implies that 

\\ag<{l + Oih^))^l^\i;{x)\^dx. (4.13) 
Again by using Schwarz's inequality, 

\ViP{x)\^dx<l [ \iV, + Vy)aix,y)\^ dxdy. (4.14) 



'•xC 



\{V, + Vy)a{x,y)f dxdy<- j _{a{. ,y)\K^'' + V{^)\a{- ,y)) dy . (4.15) 



h 

The latter expression can be bounded as 

\{\/x + \/y)a[x,y)f dxdy < 

To see this, expand a{x,y) in a Fourier series 

a{x,y)= e'P<''+y^/^ap{x - y) . (4.16) 

pG(27rZ)3 

Using that 5^(2;) = 5p(— x) for all p G (27rZ)^ we see that ()4.15p is equivalent to 

K^P'O + + 2 V{x/h) > hi'^p^ . (4.17) 

This holds, in fact, for ah p G M'^ since the left side is equal to 2K^'^ + ^h^p^ + 2£'b- 
By combining (|4.15p with (|4.13p . (|4.14p and (jiTSl) we see that ||VV'||2 is bounded by 
a constant times ||^/'||2- To conclude the uniform upper bound on the norm of ■0) 
it thus suffices to give a bound on the norm. To do this, we have to utilize the first 
term on the left side of Eq. (|4.8p . 

Eq. (j4.10p states that a can be decomposed as a = /iao?/' + Co) where oq is short for 
the operator So(— i/iV). The following lemma was proved in [6l Lemma 6]. It gives a 
lower bound on (Tr (aa)^)^/^, the 4-norm of a. This bound holds under appropriate 
decay and smoothness assumptions on oq which are satisfied in our case. (See the 
discussion at the beginning of Section O) 
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Lemma 1. For some < C < oo we have 
llalU > 



1/4 

h I \ip{x)\Ux I aoiQf jPyT - Ch^ 



-C/ji/^llV^llf (l + C/ji/^llV^IU)'^' , (4.18) 
where [•]+ = max{0, •} denotes the positive part. 

The fact that HV^-lb < C'llV'lb also imphes that HV'IU < C\\tp\\2 via Sobolev's 
inequaUty for functions on the torus. If we use also that HV'IU > llV'lb we conclude 
from (|TOD that ||a||4 > Ch^/^{\\^\\2 - CUWl^'^) for h small enough. In combination 
with (|4.8|) and (|4.13|) this implies that HV'lb < C. This shows that the norm of tp 
is indeed uniformly bounded. 

It follows that IICII2 < 0(/i^/^). To conclude the proof of KT\i . we need to show 
that also ||^||j^i < 0{h^/^). We can write 

C{x,y) = Co{x,y) + ^M_^ao{h-\x - y)) . (4.19) 

From the definition (j4.10p it follows easily that ||^o||_ffi < 0(/i^/^), using that — /i^V^ is 
relatively bounded with respect to K^'^ + V{h~^{ ■ —y)) + l. If we use the boundedness 
of the norm of -0 and, moreover, 

h^^ [ \ij{x) - iP{y)\^\Vao{h-\x - y))\^ dxdy 
= 4 |V'(p)l^ /" \Vao{x)\'^W{^hp-x)dx<0{h'^) (4.20) 

(since f \\7ao\'^\x\'^dx is finite), the bound on the norm of ^ readily follows. This 
completes the proof of Proposition [TJ □ 

Given Proposition [H the proof of the lower bound on the ground state energy is very 
similar to the corresponding one in in Sect. 7]. Let P be a state with £^'^^{T) < 0, 
and let ip be the function defined by the decomposition (14. ID . In order to be able to 
apply Theorems [2] and O we have to make sure that ■0 is in H^. For this purpose, we 
pick some e > with h < e < 1 and define ip^ via its Fourier coefficients 

ip<{p) = i^{p)9{eh-' - \p\) . (4.21) 

The function is thus smooth, and ||'0<IIh2 ^ Ceh^^ since ^p is bounded in H^. 

Let also '0> = — ip<- Since is bounded in H^, the L'^{C) norm of V'> is bounded 
by 0{he~^). We absorb the part \{il)>{x) + '4'>{y))'^o{h^'^ ~ v)) into ^, and write 

«(^,2/) = '^^^^^^p^Mh-\x-y))+a{x,y) (4.22) 



where 

(y{x, y) = i{x, y) + 2(27r)3/2/t2 (^ - ?/)) • (4-23) 
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Proposition [T] shows that < 0{h}^'^). From the bound ||^>||2 < 0{he~^) it thus 

fohows that ||(t||2 < 0{h^/'^e~^). We cannot conclude the same bound for the norm 
of cj, however. 

As in (fOjl . let A denote the operator A = -\{il)<{x)t{-ihV) + t{-ihV)ilj<{x)). 
The function t is given in (j3.ip . as in the previous section. Let be the corresponding 
Hamiltonian defined in (j2.4p . We can write 

+ \tvoH^{T-T^)+ f V{h~\x-y))\a{x,y)\Uxdy, (4.24) 

where Trg denotes again the sum of the trace per unit volume of the diagonal entries, 
as in (I33]). 

The terms in the first two lines on the right side of (j4.24p have already been calcu- 
lated. The first term is estimated in Theorem [21 and a bound on the second term was 
derived in Section [3] on the upper bound. Using the fact that the norm of ijj^ is 
uniformly bounded, as well as ||'i/'< ||_f/2 < Ce/h, we obtain the lower bound 

fBcs(r)>/,(£:GP(^^)_c(/i + e2)) 

+ ^Troi/A(r-rA)+ / V{h-\x-y))\a{x,y)\''dxdy. (4.25) 
2 JCxRS 

It remains to show that the terms in the last line of ()4.25p are negligible, i.e., of higher 
order than /i, for an appropriate choice of e ^ 1. We shall use the following lemma, 
whose proof is inspired by Lemma 1] . 

Lemma 2. For all Q <T < 1 with (— + 1)7 trace class, we have 

Tro Ha(T - La) > Tr (r - La) I^^aI (r - La) • (4.26) 

Proof. Recall that Fa is the projection onto the negative spectral subspace of H/\. 
Moreover, for any operator ^, Tro^ = TrFo^iro + Tr (1 - Fo)^(l - Fq). A simple 
calculation shows that 

ToHAiT - Fa)Fo + (1 - Fo)i^A(r - Fa)(1 - Fq) 
= H+r{l - Fa) +H^{1- F)Fa - Ei - E2 (4.27) 

with denoting the positive and negative parts of H/^, respectively, 

^1 = (ro-FA)i^A(r-FA)(l-2Fo) (4.28) 

and 

i?2 = |^A|(r-FA)(Fo-FA). (4.29) 

It is easy to see that (Fo-Fa)IFoI^/^ and |/7o|^/^(r - Fa) are Hilbert-Schmidt; since 
A is bounded, also (Fq - FA)|i/A|^/^ and |i/A|^/^(r - Fa) are Hilbert-Schmidt. Hence 
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El is trace class and, by cyclicity, its trace is equal to the one of 

El = v^(r - rA)(i - 2ro)(ro - rA)^^ 
- y/H^{T - Ta){i - 2ro)(ro - rA)^^ 
= -/^(r - rA)(ro - Ta)/^ - /^(r - rA)(ro - rA)/^^ • (4-30) 

Via the Floquet decomposition, can be written as a direct integral of opera- 
tors on LP'{C) each of which has discrete spectrum. Since we know, a priori, that 
(j4.27p is trace class, we can evaluate the trace in the basis given by -f^A- With this 
understanding of the trace, we have Tr [Ei + £"2] = 0, and thus 

Tro H^{T - Ta) = Tr [H+T{1 - Ta) + H^{1 - T)T^] 



Tr 



'hXtJhX + Jhi{i-t)Jhi 



(4.31) 



The operators on the last line are positive, hence they are trace class. Estimating 
r > and 1 — r > (1 — r)^, respectively, gives the desired bound 



Troi/A(r-rA) >Tr 
= Tr 



(4.32) 
□ 



which agrees with the right side of (j4.26p . 

An application of Schwarz's inequality yields 

Hi>il-r^)Hl-r]~'\\A\\l (4.33) 

for any 7] > 0. Schwarz's inequality can also be used to obtain a lower bound on Hq. 
For any < S < 1, 



-ihV + hA{x) f + £b + h^W{x)\ > (1 - 6f [-/i^V^ + £b 



(-i/i^V • A{x) - ih'^A{x) ■ + {W{x) + ^(5 



^2^2■ 



(4.34) 



We can further bound 

{-ih^V ■ A{x) - ih^A{x) ■ V)"^ 

= (-2i/i2 V • A{x) + ih^ div A{x)) {-2ih'^A{x) -V -ih^ div A{x)) 
< 8/i^V • A{x)A{x) • V + 2/i^ {dwA{x)f . 



(4.35) 



Since A is by assumption, this is bounded from above by C/i^(— + 1). Choosing 
6 = 0(h), we thus conclude that 

> {I - 0{h))[-h''V^ + Ei,f 



The operator monotonicity of the square root implies that 



(4.36) 
(4.37) 
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Using again (ICTI) and the fact that > (A^'^)^ > 0(1), the choice rj = 0(||A||oo) 
gives 

\Ha\ > (1 - 0{h + ||A||oo))-?ir°'° IIc2 . (4.38) 
In particular, we infer from (j4.26p that 

^TtHa{T-Ta) > (l-O(/i + ||A||oo))Tri^°'0(a-aA)(a-aA), (4.39) 

where a a denotes again the upper off-diagonal entry of Ta- From the definition of A, 
we see that 

||A||oo < /i||V'<l|oo|N|oo- (4.40) 
Moreover, since the Fourier transform of ip^ is supported in the ball \p\ < e/h, 

||V'<IIoo<^|^<(p)|<||V<IIhi(c) f E TT^] <cV^, (4.41) 

and hence ||A||oo < 0{e^/'^h^/^). 

Recall the decomposition ()4.22p of a. We decompose oa in a similar way, and define 
4> by 

OA = I {iJ<ix)ao{-ihV) + ao{-ihV)i;<{x)) + (t> . (4.42) 
In particular, we have 

a — OA = a — 4> ■ (4.43) 
Since ||^<||j:/2 < 0{e/h), Theorem [3] implies that ||</>||/i-i < 0{eh^/'^). From the posi- 
tivity of K^'^ we conclude that 

TV K^'° (o- - 0) (a - (^) > Tr K^'^aa - 2 Re Tr K^'^a^^ . (4.44) 

The terms quadratic in a are thus 

(1 -(5)Trir°'VCT+ / V{h-^{x-y))\a{x,y)\'^dxdy (4.45) 
JcxR^ 

with 5 = 0{h + ||A||oo) = 0(ei/2/iV2). pi^k some 5>0 with 5 + 6 < 1/2, and write 
(1 - 6)K^'° + V = 5K°'^ +(l- 26 - 26^ {K^'° + V) + [6 + 6^ (i^°'° + 2V) 

> -c(6 + 6) , (4.46) 



where we have used that V is relatively form-bounded with respect to K^'^ to bound 
the last term. Hence ()4.45p is bounded from below by 

'^{hWm -{C-E^ + l)\\a\\l) - C5Ml . (4.47) 

Recall that ||cr||2 < Oih}!"^ j e). We shall choose = if the first parenthesis on the 
right side of (|4.47p is less than ^HcH^i (and, in particular, if it is negative), while 
b = 0(1) in the opposite case, i.e., when ||cr||^^ > 2{C — E^, + l)||cj||2. In the latter 
case we shall have the positive term (5||(t||^i/2 at our disposal, which will be used in 
([130]) below. 

We are left with estimating the last term in (j4.44p . which is linear in a. Recall from 
(j4.23p that o" is a sum of two terms, and cr — ^, where the latter is proportional to 
^/'>, and ll^lliji < 0{h^/'^) independently of e. Moreover, as the proof of Theorem [3] 
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shows, (j) is the sum of two terms, r^i and (p — r]i, with rji defined in ()2.14p (with ip 
replaced by '0<) and ||</' — ?/i H/fi < 0(/i^/^). Now 

Tr Er°'° (o- - ^) f/i = (4.48) 

as can be seen by writing out the trace in momentum space and using that and 
■0> have disjoint support. Hence 

Re Tr^O'Vc^ <Cm\H4<P\\m + ||ct||^i||0 - r/iH^i) 

< 0(e/i) + 0(/i3/2)||a||^i . (4.49) 

In the case < C||(t||2 (corresponding to (5 = above) we can further bound 

llfll/i-i < 0(/i^/^/e). In the opposite case, where 5 = 0(1), we can use the positive 
term 5||cr||^i/2 from before and bound 

^-\\a\\l,-Oih'/')\\a\\Hr>-0{h'), (4.50) 

which thus leads to an even better bound. 

In combination with ()4.25p these bounds show that 

£^^^(r) > h {£^^{^<) - Ce) (4.51) 

where 

h h^/'^ 

e = /i + + £+- + — -. (4.52) 

The choice e = h^/^ leads to e < Ch^/^. 

The completes the lower bound to the BCS energy. The statement (jl.9p about 
approximate minimizers follows immediately from ()4.5ip and (j4.22p . 
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